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0> : ABSTRACT 

p 

" We use near-infrared spectroscopic data from the inner few hundred parsecs of 

^SJ . a sample of 47 active galaxies to investigate possible correlations between the stellar 

velocity dispersion (cr^), obtained from the fit of the K-band CO stellar absorption 
. bands, and the gas velocity dispersion (cr) obtained from the fit of the emission- 

5_j ; line profiles of [S in]/l0.953//m, [Fen]/11.257;um, [Fe ii]/11.644/zm and H2i2.122//m. 

While no correlations with cr^ were found for H2 and [S m], a good correlation was 
found for the two [Fell] emission lines, expressed by the linear fit cr,^ - 95.4 + 
16.1 + (0.25 + 0.08) X a-[psii]. Excluding barred objects from the sample a better 
correlation is found between 0-^, and cr[pe n] , with a correlation coefficient of = 0.80 
and fitted by the following relation: cr,^ - 57.9 + 23.5 + (0.42 + 0.10) x o-fpen]. 
This correlation can be used to estimate 0-^, in cases it cannot be directly measured 
and the [Fe 11] emission lines are present in the spectra, allowing to obtain the mass 
of the supermassive black hole (SMBH) from the M, - o-^, relation. The scatter 
from a one-to-one relationship between cTi, and its value derived from cr[pe nj using 
the equation above for our sample is 0.07 dex, which is smaller than that obtained 
in previous studies which use cr[oiu] in the optical as a proxy for (t*. The use of 
o"[Feii] in the near-IR instead of crjoiu] in the optical is a valuable option for cases in 
which optical spectra are not available or are obscured, as is the case of many AGN. 
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The comparison between the SMBH masses obtained using the M. - cr^ relation in 
which (Ti, was directly measured with those derived from cr^pe nj reveals only a small 
average difference of AlogM, = 0.02 with a scatter of 0.32 dex for the complete 
sample and AlogM, - 0.00 with a scatter of 0.28 dex for a sub-sample excluding 
barred galaxies. 
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1 INTRODUCTION 



In the present paradigm of galaxy evolution, most galaxies which form a bulge also forrn a su 



permassive black hole (SMBH) in their nuclei (e.g. 



Ferrarese & Merrit 



2000; 



Gebhardt et al. 



Magorrian et al.l 



1998 



Richstone et al. 



1998 



2000|) . The central SMBH seems to play a fundamental 
role in the galaxy evolution and cosmological simulations without considering the presence of 
a SMBH and its associated feedback predict masses for the galaxies much higher than those ob- 



serve d (jpi Matteo. Springel & Hernquist 2005 ; Springel. Pi Matteo & Hemquist 2005 



Bower et al. 



20061). In a scenario of co-evolution of the SMBH and its host galaxy, mass accretion to the central 
region of the galaxy leads to the growth of the galaxy bulge, while the feeding of the SMBH trig- 
gers episodes of nuclear activity which results in feedback in the form of radiation pressure and 
mass ejections from the accretion disk surrounding the SMBH. This episo dic feedback rnay hal t 



the mass accretion to the galaxy preventing its growth in the active phase (|Nemmen et al 



2007h 



This co-evolution may be the mechanism which leads to the empirical relation between the mas s 



of the SMB H and the stellar velocity dispersion of the bulge M, - ct-i,, [[Ferrarese & Ford 



but see also 



Jahnke & Maccio 



most) : 



mm]- 

The M,-(T-t, relation has been extensively used to estimate the mass of SMBHs from the stellar 
kinematics, as direct determinations for the SMBH masses only are poss ible for the closest galax - 



Ferrarese & Ford 



20051) 



ies for which the radius of influence of the SMBH can be resolved (e.g 
Although allowing to estimate the masses of the SMBHs for a large number of galaxies, the use of 
the M, - cr^ relation requires the measurement of cr^, which is not always easy to obtain, particu- 
larly in active galaxies, where the AGN continuum dilutes the stellar absorption lines. In order to 



overcome this difliculty, a number of scaling relation s using the width and 



sion 


ines to determine M, have been proposed (e.g. 


Boot 


1 & Schave 


2011 


; > 


^u 


20( 


)9; 


Peterson 


2008 


Salviander et al. 


20061; 


Vestergaard & Peterson 


2006; 


Greene & hJI2006 


,2005; 


Kaspi et al. 



uminosities of emis- 
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2005 



Onken et al. 



2004 



Nelson & Whittle 



19961) . Nevertheless, most of these relations are for the 



optical domain of the electromagnetic spectrum. With the improvement of infrared (IR) detectors, 
IR spectra of many AGNs have become recently available, and have the advantage of being less 
affected by reddening than optical spectra. In the present paper, we investigate the possibility of 
using the widths of emission lines in the near-IR as proxies for cr^,. 

Recent studies by our group, using near-infrared (hereafter near-IR) integral field spectroscopy 
of active galaxies, have allowed the mapping of the flux distributions and kinematics of the molec- 
ular (H2) and ionized gas. We have found, in particular, that the H2 usually shows small velocity 
dispersions and a velocity field dominated by rot ation, while the ionized gas shows higher velocity 



disper sions and is not dominated by rotation (e.g. 



Riffel et al. 



2008 . 



20091 : 



Riffel & Storchi-Bergmann 



2011alJbh . The kinematics and flux distributions are also consistent with a location for the H2 gas in 



the galaxy plane, while the ionized gas, and, in particular [Fe 11] extends to high galactic latitudes. 

In the present paper we investigate the correlation between the gas and stellar kinematics de- 
rived from near-IR spectroscopy, with the goal of looking for a "proxy" for cr^ among the brightest 
emission lines in this wavelength range. This paper is organized as follows: in Section 2, we de- 
scribe the sample and the observational data; in Section 3 we describe the methods used to measure 
the stellar and gaseous velocity dispersion. The results are presented in Section 4 and discussed in 
Section 5, while the conclusions are shown in Section 6. 



2 OBSERVATIONAL DATA 



The spectroscopic data used in this work are f rom 



Rodriguez- Ardila. Contini & ViegasI (|2005|) and [Rodriguez- Ardila et al 



Riffel. Rodriguez-Ardila & Pastorizal (120061) : 



(120041) . The sample com- 



prises 47 active galaxies with a range of activity types, and the spectra cover, on average, the inner 
300 pc radius of the galaxies. 

The near-infrared spectra were obtained with the NASA 3 m Infrared Telescope Facility (IRTF), 
using the SpeX spectrograph in the short cross-dispersed mode (SXD, 0.8-2.4/im). The detec- 
tor employed was a 1024x1024 ALADDIN 3 InSb array with a spatial scale of 0.15"/pixel. A 
0.8"xl5" slit was used and the spectral resolution is 300kms"\ obtained from the measurement of 
the full width at half maximum (FWHM) of Arc lamp lines, or 127 km s"^ in cr. The data reduction 
followed standard procedures. For more details on the instrument al configuration, data reduction. 



calibra tion processes and details of the extraction of the spectra see lRiffel. Rodriguez-Ardila & Pastoriza 



(2006). 
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The above sample was chosen for this work because it is an unique dataset of near-IR spec- 
troscopy of active galaxies, observed with the same instrumental setup (thus with no instrumental 
bias), covering the near-IR Z, J, H and K bands, and including several emission and absorption 
features allowing the comparison of the stellar and gas kinematics. 



3 METHODS 



In order to obtain the gaseous velocity dispersion cr, we fitted the emission-line profiles of [S m]/l0.953/im, 
[Feii]/ll.257yum, [Feii]/ll.644//m and H2/l2.122/im by single Gaussian curves and adopted as the 

measured velocity dispersion the cr of the Gaussian. These emission l ines have been chosen be- 

cause they are the strongest in the near-IR spectra of active galaxies (e.g. lRiff'el. Rodriguez- Ardila & Pastoriza 
20061) . We excluded the H and He recombination lines due to uncertainties in the fit for type 1 
objects, for which is not always easy to separate the narrow from the broad components. The fit- 



ting of the emission-line profiles was done by adapting the profit routine (IRiflrel 



2010h . which 



outputs the emission-line flux, the centroid velocity, the velocity dispersion and the uncertain- 
ties for each of these parameters. The velocity dispersion was then corrected for the instrumental 
(Tinst =127 kms \ which was subtracted in quadrature from the cr obtained from the fit of the 
Gaussians to the line profiles. 

We me jasured the stellar velocity dis version (cr^) using the penalized Pixel-Fitting (pPXF) 
method of Capp ellari & EmsellemI (|2004() in order to fit the CO absorptions bands at ~ 2.3/im 
in the K-band. The pPXF method requires the use of stellar spectra as templates. We used for 
this th e spectra of 60 late-type stars, 40 of the m from the Gemini Near-IR Late-type stellar li- 
brary (jWinge. Riff"el & Storchi-Bergmannll2009h and the remaining 20 spectra are from stars with 
public NIFS observations in the Gemini data archive (Diniz et al., 2012 in preparation). The un- 
certainties on the measurements o f cr^ were estimated using 100 Monte-Carlo iterations as in 
Riffel & Storchi-BergmannI (1201 Ibl) . 

In order to illustrate our procedures, we show in Fig. [U sample fits of the CO absorption band 
heads at 2.3 jjm using the ppxf method, as well as fits of the emission-line profiles using profit for 
the spectrum of the galaxy NGC5929. 



4 RESULTS 

The resulting measurements for the stellar and gas velocity dispersions for the galaxies of our 
sample are shown in Table [TJ The dashes in the Table are due to the fact that for a few objects we 
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Figure 2. Compaiison between the gas velocity dispersions obtained from the H2/12.122/im emission line (ctht) and the stellar velocity dispersions 
from the CO stellar absorptions at ~2.3,um (o"^). Points with no error bars in one or both axes represent measurements that are unresolved by our 
observations and should be considered as upper limits. 

were not able to measure one or more values due to the absence of the absorption/emission lines 
or due to a low signal-to-noise ratio. 

We have looked for correlations between the stellar and the gaseous velocity dispersions using 
the values of Table [T] to construct the graphs of Figures [31 [2l and |4l We have used the IDL routine 
R CORRELATE to obtain the Spearman correlation coefficient R for each graph. Figure [2] shows cr^f vs. 
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Table 1. Stellar and gas velocity dispersions for the galaxies of the sample. 



Object 


Hubble Type'' 


Nuclear Activity" 


o-y, (km s ' ) 


(TH2 (km s ') 


o-[Feii] (km s ^) 


o"[Siii] (kms ') 


iVilvJV J 


Sbc 


syi 


<127 


<127 


180+9 


278+14 


NGC34 


so/a 


oyz. 


160+23 


198+17 


227+19 


215±18 




so/a 


Sy2 




<127 


263+17 


245+16 


MRK'QQ^ 

iVil\.J\.?7 7-J 


Od 




132+14 




214+12 


287+16 




SBO/a 


Sv9 
ayz 


130±21 


<127 


274+17 


229+14 




SBO 


Sv9 
syz 


<127 


<127 


172+12 


211+15 




SB(s)b 


Syl 


165+13 


<127 


166±12 


381+28 


Nnn 144. 




Sy2 


7(16+40 


175+13 


230+18 


233+18 


MRKIOfifi 


so/a 


Sy2 


<127 


<127 


206±15 


234+17 


i^l vJv- i Z. / J 


SO 


Sy2 




173+21 


313+38 


564+68 




SB(s)c 


Sb 


<127 


<127 


203+12 


217+13 




SBO/a 


Syl 


162+11 


191+22 


162+19 


201+23 




E-SO 


Sy2 


184+16 


<127 


248+17 


305+21 




SO 


oyz< 


<127 


161+1 1 


233+18 


253+18 


MR K 1 9 1 


S? 


Sy2 


181+24 


<127 


316+32 


342+35 


MR If 194. 


S? 


NLSl 


222+23 


235+27 


335+39 


244+29 


MRK19'^0 


E-SO 


NLSl 






<127 


370+35 


"MnP^997 


SB(s)a 


Syl 


128+3 


157+24 


314+48 


275+42 


HI 143-1Q2 

±111 '-x^ 1 y ^ 




Syl 






17Q+28 


410+68 




OlJ(^l _/uu 


Sb 


142+18 


<127 


130+9 


153+10 


pn 1 1 9f^ 041 


c 


QSO 




<127 


254+34 






OlJl^l f» )\J\^ 


NLSl 


<127 


<127 


142+6 


215+9 




vjJJl^l a fa\J 


Syl 


136+19 


<127 


194+41 


231+48 


MRK7fSfi 

IVIIVIX. 1 \J\J 


SB(s)a 


NLSl 


<127 


<127 


133+20 


188±19 


NGC4748 


S? 


NLSl 


<127 


170±15 


171+21 


273±23 


TONSO! 56 

1 vyi 1 ljU 1 K-'U 




QSO 






<127 


712+200 


MRK27Q 

iviiviv^ / y 


SO 


NLSl 


138+12 


<127 


209±49 


335+79 




so/a 


Syl 


<127 


<127 


165+43 


213+55 


MRK47R 

J.V11V1VT / o 


Sc 


NLSl 




161+22 


171+23 


390+53 


i'IVJ>_J / Z.O 


SBa 


oyz< 


<127 


132+23 


<127 


223±38 










< 127 






MRK6R4 


Sab 


Syl 


170+17 


351+18 




540+27 


MRK'901 


SBa 


NLSl 




129+14 


144+16 


167+18 


MRK49'? 


SBb 


NLSl 




<127 


232±12 


412+21 


NGC5929 


Sa 


Sy2 


158±23 


<127 


208±26 


181±23 


NGC5953 


so/a 


Sy2 


149±4 


187±8 


289±13 


289±13 


PG1612+261 




QSO 


166±18 




266±53 


262±52 


MRK504 


S? 


NLSl 


<127 


<127 




333±38 


3C351 




BLRG 








187±14 


ARP102B 


EO 


Syl 


<127 


155±9 


195±10 


299±16 


IH 1934-063 




NLSl 


184±9 


191±8 


238±11 


291±13 


MRK509 


E-S? 


Syl 


<127 




159±39 


370±91 


1H2107-097 




Syl 


<127 




<127 


256±56 


ARK564 


SBc 


NLSl 




180±10 


176±10 


260±15 


NGC7469 


SBbc 


Syl 






161±15 


283±18 


NGC7682 


SBab 


Sy2 


183±40 


<127 


198±27 


205±28 


NGC7714 


SB(s)b 


HII 


<127 


<127 


165±13 


184±14 



" The Hubble type and Nuclear Activity were taken NASA/IPAC Extragalactic Database (NED) and Hyperleda Database ^Paturd^^al^^MTO). 

cth, . The range of the cr^ and values is approximately the same, something we have noticed in 
our previous stu dies of individual galaxi e s using integral field spectroscopy of the inner hundreds 



of parsecs (e.g. 



Rifi^el et al. 



2008 



20091 : 



Riffel & Storchi-Bergmann 



2011a 



b|). Nevertheless, we 



have obtained only a very weak correlation between cth^ and cr^, with R = 0.35, but we note that 
the H2 line is unresolved for several objects. 

The relation between (Ti, and cr[pen] is presented in Figure |3l showing that cr[pen] is usually 
higher than cr^, which is also in agreement with the results from the integral field spectroscopic 
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Figure 3. As Fig.|2]but for gas velocity dispersion derived using [Fen]/il.257/jm. Tlie daslied line represents tlie best linear fit of the data, given 
by: a■^, = 95.4 ± 16.1 + (0.25 ± 0.08) x o-[Feii]. 
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Figure 4. As Fig.|2]but for gas velocity dispersion derived using [S ni]/l0.953^m. 



Studies above. A better correlation is observed between cr^f and crpeii] than with cth,, corresponding 
to a correlation coefficient R = 0.56 obtained as described above. We fitted the data by a linear 
equation of the form cr* = A + Bx cr [Peii] using the IDL routine linmix_err, which uses a Bayesian 
approach to linear r egression wi th errors in both variables and takes into account upper limits for 



the measurementes ( Kelly 



20071). The best fit to the data is given by 
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200 300 

<^rFe lllA1.25«m(kn^/s) 

Figure 5. cr[Fcn],i 1.644 vs. o"[Peii],u.257- The solid line represents a one-to-one relationship. The top panel shows the difference between the values 
of o"[peii]^i 644 and 0"[p5iij^i 257 defined as D(<t) = (T[peii]^i 541 - (TiVcllui.lSiy as well as the average value of this difference (D). 



(T^ = 95.4 + 16.1 + (0.25 ± 0.08) x (^[Feii] (1) 

which is shown as a dashed line in Fig. [31 

Finally, the relation between cr^^ and cr[sni] is shown in Figure HI resulting in a correlation 
coefficient R = 0.32, suggesting only a very weak correlation. This figure also shows that cr[siii] is 
larger than cr* by more then a hundred km s~\ on average. 



5 DISCUSSION 



The use of the velocity dispersion from the Narrow Line Region emission lines as a proxy for 
in order to obtain an estimate for the SMBH mass via the Mbh -c* relation in active galaxies is not 
new. It has been previou sly used in the optical domain, where the emission line most cornmon 



used is [Oiii]i5007 (e.g. 



Wu 



2009 



Salviander et al. 



20061 : 



Kaspi et al. 



20051 : 



Onken et al. 



2004) 



This emission line has been used instead of cr^ because in active galaxies cr^ cannot be easily 
measured due to dilution of the stellar absorption lines by the AGN continuum or its scattered 
light. 

In this paper we present an alternative to be used in the near-IR. As shown above, we found 
a correlation between cr^ and cr[Feii], indicating that the latter can be used to estimate cr^ using 
equation [B [Fe ii] has two similarly strong emission lines which can be observed in the near-IR: 
[Feii]/1 1.257 in the J band and [Feii]/l 1.644 in the H band. In Figure |5] we present a comparison 
between the cr of these two lines, where the solid line shows an one-to-one relationship. This 
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comparison shows that the width of these lines is the same within the errors, with a mean difference 
of cr[Feii]ii.644-C[Feii]ii.257 = 8 km s"' and a scatter of 47 km as seen in the top panel of Fig. [51 
This scatter may be partially due to the fact that [Feii]/I1.644 is close in wavelength to Brackett 12, 
which usually appears in absorption and may affect the measurement of [Fe ii] line. 

Why is CTffg//] better correlated with cr^, th an cth^ ? As pointed out in the Introduction, our previ 



ous stu dies using Integral field spectroscopy (IRiffel et al. 



2008 



20091 : 



Riffel & Storchi-Bergmann 



2011alJbh showed that the H2 kinematics frequently shows a rotation pattern and a smaller velocity 



dispersion than that of the ionized gas. This also led to the conclusion that the H2 gas was more 
restricted to the galaxy plane, while the ionized gas - and in particular [Fe 11] - extended to higher 
galactic latitudes. The integrated value of cr^ from the nuclear region of galaxies is dominated by 
the contribution of bulge stars, which are not restricted to the plane, showing a "hotter" kinemat- 
ics. Thus, it can be understood that the velocity dispersion of gas which is restricted to the plane 
does not correlate with that of bulge stars, while that of gas extending to higher latitudes, similar 
to those of the bulge - such as the [Fe 11] emitting gas - is correlated to that of the bulge stars. The 
higher values of cr[Feii] relative to cr^^ are probably due to extra heating provided by a nuclear AGN 
outflow. 

Our results can be compared with previous ones in the optical usin g the cr of the [O m 1 /i5007 



Nelson & Whittle 



1996) 



Onken et al. 



emission line as a proxy for cr^. For a sample of 66 Seyfert galaxies 
found a scatter of 0.20 dex around a one-to-one relation between cr^, and cr^oni], while 
([2004 ) found a smaller scatter of 0. 1 5 dex for a sample of 1 6 AGNs, for which 25% of their sources 
have CTfoiii] deviating by more than 0.20 dex from the values expected based on their cr^. We found 
a scatter of 0.07 dex between the values obtained via equation [T] and the measured values for cr^, 
which is thus smaller than that for crjoin]. 

A cautionary note is the observation of recent spatially resolved studies (e.g. our previous stud- 
ies already mentioned) that the [Feu] emission originates at least in part in outflowing gas. Thus, 
the width of the line is not only due to orbital motion in the gravitational potential of the galaxy, 
but is also due to broadening by the outflow, what is consistent with the observation that o-^p^u] 



is larger than cr^^. Similar outflovys - most probably the same - are observed in the [O m 



ting gas (e.g. 


Fischer et al. 


2011 


Das. Crenshaw & Kraemei 


2007 



2010; 



Crenshaw et al. 



Das et al 



2010 



2009 



Crenshaw & Kraemer 



emit- 



2007 



20051). Nevertheless, this line has been frequently used 



as a proxy for cr^ as discussed above, due to the lack of a better indicator. Our argument is that the 
cr [Fell] is at least as good cr* proxy as cr^oiii], and can be used when the latter is not available. 
We thus propose the use of cr[Feii] to obtain cr* via eq.[I]in cases for which it is not possible to 
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Figure 6. As Fig.[2]but for gas velocity dispersion derived using [Fe ii]/t 1 .251 lira for tlie sub-sample of unbarred galaxies. The dashed line represents 
the best linear fit of the data, given by: cr* = 57.9 ± 23.5 + (0.42 ± 0.10) X (T[Feii]. 

measure the stellar kinematics of the galaxy, and the optical spectrum is obscured or not available, 
so that the [O III]/15007 emission line is also not available. Nevertheless, this suggestion should be 
used with care, since the M, - cr^ is calibrated from a parent sample of mostly early-type galaxies 
and, as seen in Table[i| most of the objects of our sample are late-type. Late-type galaxies can 
have distinct cr^ values than early-type galaxies, since the orbits of the stars in a disk are different 
than the orbits of the stars in the bulge. Additionally, some of the galaxies of our sample have bars, 
circumnuclear star forming rings or nuclear starbursts or even are classified as peculiar objects and 
thus the cr^ measured for these objects could be very difl'erent than those for the classical bulge, 
used to calibrate the M,* - cji, relationship. 

5.1 The effect of galaxy morphology on the M, - cr relation 

As pointed out above, the M, -cr^ relation is calibrated using a parent sample of mostly early-type 
galaxies. But Table[T] shows that about 30% of the galaxies of our sample are late-type, which can 
have distinct cr^ values from those of early-type galaxies, since the orbits of the stars in a disk 
(which dominate in late-type galaxies) are distinct from those in a bulge (which dominate in early- 
type galaxies). Additionally, « 30% of the galaxies of our sample have bars, and another 30% have 
uncertain classifications and are peculiar objects which may not obey the the M, - cr^, relationship. 
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Xiad (|2011|) investigated the M, - cr^ relation for late-type galaxies for a sample of 93 objects 
with Seyfet 1 nucleus. They examined the M, - cr* relationship for subsamples of barred and 
unbarred host galaxies and found no difference in slope. They only found a mild offset in the 
relation between low- and high-inclination disk galaxies, with the latter tending to have larger cr^ 
for a given value of the black hole mass. 

The M, - cr^ re lationship for galaxies of different Hubble types have also been studied by 



Graham et al. 



(|2011l) using a sample of 64 galaxies. They found that restricting the sample only 
to elliptical galaxies, or only to non-barred galaxies result in tighter relations (with less scatter) 
and a smaller slope than when using the full sample of galaxies. The M, - cr^ relation obtained 
when the sample is restricted to barred galaxies only lies « 0.45 dex bellow the relation obtained 
for elliptical and non-barred galaxies. 

In order to investigate the effect of the presence of a bar in the M, - cr[Feii] relation of Fig. [31 
we divided our sample into two sub-samples: one composed of barred galaxies only and the other 
of unbarred galaxies. We found a much better correlation between crpeii] and cr^^ for the unbarred 
galaxies than for the total sample, as illustrated in Figj6l The correlation coefficient is R = 0.80 
and a linear regression to the relation is given by 

cr^ = 57.9 + 23.5 + (0.42 ± 0.10) x (r[Feii]. (2) 

On the other hand, no correlation was found for the barred galaxies, for which the correlation 
coefficient between o-^p^u^ and cr^ is only R = 0.20. 

Approximately half of our sample is composed by early-type galaxies (half of which are 
barred), the remainder being of late-type galaxies and galaxies with uncertain morphology (usu- 
ally because they are distant and compact). But the number of galaxies of our sample is not large 
enough to further investigate the effect of morphology (besides the effect of bars discussed above) 
on the (Ti, - CTfpeii] relation. This will be possible only when more near-IR AGN spectra become 
available, and the equations [Hand [2] could be better calibrated using large non-biased samples of 
galaxies. 

On the other hand, when applying the cr^ - o-^peu] relation to distant galaxies, it will be hard to 
ascertain a Hubble type to such galaxies, and it will probably better to just use the relation for the 
whole sample (equation [Tj), in spite of the fact that a relation for restricted Hubble types (e.g., for 
galaxies without bars) shows less scatter. 

Many present/future missions are discovering/will discover large numbers of obscured AGN, 
such as the Wise mission 
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\http : I lwww.nasa.gov Imissiori-pages IWIS E Inews I wise2Q\2Q'&29.htmJ) or the Vista Variables 
in the Via Lactea survey 

\http : / /mwm.astro.puc.cl/mw/index.ph p/MainJ'age), which are potential samples to benefit 
from the relation we have found in order to obtain the SMBH masses, if [Fe ii] emission lines are 
present in the spectra. 



5.2 Estimating SMBH masses using crpen] 

We now evaluate the effect of the scatter introduced by the use of eq.[T]in the value of SMBH 
masses obtained via the M, - cr* relation, which result in higher uncertainties in M, than those 
obtain ed by using cr* directly. In order to do this, we use the relation below from iGraham et al. 

honk 

log(M./Mo) = (8.13 + 0.05) + (5.13 + 0.34)log[cr^/200kms-i]. (3) 

to obtain M,/Mq, using for cr^ first its measurement from the stellar kinematics and then the 
value derived from cr[pgn] using eq.\T\ We compare the two values in Figure |7] (log M(cr^) vs. 
log M(cr[pen])), which shows a good agreement between them, with an average difference of AlogM, = 
logM(cr[pen)-logM(cr+) = 0.02+0.32. In the same figure, we show also the masses for the SMBH 
obtained for the sub-sample of unbarred galaxies as open symbols, using the eq.[2]to obtain cr*. 
There is no diference in the average value of the SMBH masses obtained, and the mean scatter is 
0.28 dex, which is somewhat smaller than the one for the complete sample. 



6 CONCLUSIONS 

We have used near-infrared spectroscopic data of a sample of 47 active galaxies in order to inves- 
tigate possible correlations between the stellar velocity dispersion (cr^^), obtained from the fit of 
the K-band CO absorption band-heads, and the gas velocity dispersion (cr) obtained from the fit 
of the profiles of the [S iii]i0.95332jum, [FeII]/ll.25702jum, [Fen]il.644//m and H2/l2.12182//m 
emission lines. The main conclusions of the present paper are: 

• Very weak correlations are found between cr^^ and both 0-^2 and cr[sni]; 

• The best correlation is found for the [Fell] emitting gas with R = 0.58 for the Spearman 
rank correlation coefficient between cr* and cr[Feii] (both the /ll.257/zm and /ll.644jum emission 
lines can be used). A better correlation is found if we exclude the barred galaxies from the sample 
(R=0.80), while no correlation is found for the sub-sample of barred galaxies. 
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D = 0.02 a = 0.32 /unbarred D = 0.00 cj=0.2B 




logW^rFelll)/^©] 

Figure 7. Comparison between SMBH mass values obtained directly from the relation with cr* (y-axis) and using a-yp^w] (x-axis) to obtain rr*. 
The filled squares are for the complete sample, for which the eq.[T]was used and the opened squares points of for unbarred galaxies, for which the 
eq.|2]was used. The error bars shown include the uncertainties in the calibration of the M. - it-i, relation. 



• cr[Feii] can thus be used to estimate cr* for objects for which the stellar velocity dispersion 
cannot be measured or is unknown. The best fit of the data is given by the equation cTi, = 95.4 + 
16.1 + (0.25 + 0.08) X cTfFeii] for the complete sample and cr^ = 57.9 + 23.5 + (0.42 + 0.10) x o-[Feii] 
for the sub- sample of unbarred galaxies. 

• The equations above should be improved and re-calibrated when larger and non-biased sam- 
ples of near-IR spectra of AGN become available. 

• The scatter from a one-to-one relationship between cr^ and its value derived from crfpen] 
using the equation above for our sample is 0.07 dex, which is smaller than the scatter of previous 
relations using CTfom] in the optical as a proxy for cr^. 

• The use of cr[Feii] in the near-IR instead of c7-[oiii] in the optical is particularly important for 
cases in which the optical spectra is not available or is obscured, as is the case of many AGN. 

• The comparison of the masses for SMBHs obtained from the direct use of cr^ in the M, - cr^ 
relation with those using (T[Feii] to obtain cr^ reveals only a small average difference of AlogM, = 
0.02 + 0.32 for the complete sample and AlogM, = 0.00 ± 0.28 excluding barred galaxies from the 
sample. 
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